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We investigate the fermion creation in quantum kinetic theory by ap- 
plying "oscillator representation" approach, which was earlier developed for 
bosonic systems. We show that in some particular cases (Yukawa-like inter- 
action, fixed direction of external vector field) resulting Kinetic Equation 
(KE) reduces to KE obtained by time-dependent Bogoliubov transforma- 
tion method. We conclude "oscillator representation" approach to be more 
universal for the derivation of quantum transport equations in strong space- 
homogeneous time-dependent fields. We discuss some possible applications 
of obtained KE to cosmology and particle production in strong laser fields. 

PACS numbers: ll.10.Ef, 25.75.Dw 

1. Introduction 

Spontaneous particle production of pairs under the action of strong 
external fields attracts a lot of attention since it is a multi-particle non- 
perturbative effect which is still lacking experimental verification. Accord- 
ing to the famous Schwinger formula [T] the probability of particle creation 
in the constant electromagnetic field becomes essential if the field reaches 
the critical value E c = m 2 /e. Considering electron-positron pair creation 
in laser field we may conclude that one must have unattainable static field 
E^ +e _ = 1.32 • 10 16 V/cm. However theoretical and numerical results ob- 
tained in recent papers (e.g. EJ HI EJ E]) show that the probability of 
particle creation in time- dependent fields can be gradually enhanced in con- 
trast to static ones. Also recent developments in laser technology (method 
of chirped pulse amplification [7] ) and underway construction of X-ray free 
electron lasers will allow an experimental verification of spontaneous particle 
creation in the nearest future. 

Fermion production in external time-dependent fields was considered 
in a number of papers [U [TO] lllj . The method of the time-dependent 
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Bogoliubov transformation used there faces difficulties if the external field 
has more than one non-zero component, thus KE was obtained only for the 
linear polarized (also called flux-tube geometry) field (A^ = (0, 0, 0, A(t) = 
^3 (^)))- A s was shown earlier for boson creation [15] in periodic laser field, 
pair production is efficient if the field has rather the circular polarisation 
than the linear on^|. 

The approach presented in current paper resolves the problem of fermion 
creation in an arbitrary polarised field (e.g. in Hamilton gauge A^ = 
(0, Ai(t),A2(t),A^(t)) ). This result can be of great importance for the 
theoretical predictions about dilepton and photon yield in high-intensity 
laser experiments [131 14] . 

As already mentioned, the classical approach to the description of the 
particle creation in the external fields or in curved space is based on the 
time-dependent Bogoliubov transformation [8] [9]. Concentrating on fermion 
fields we will briefly summarise this approach. First of all let us suppose 
that in the asymptotic states t — ► ±oo the external field vanishes. It means 
that for the limit t — > ±oo the interpretation in terms of particles and 
antiparticles becomes possible. Thus we can introduce the orthogonal set 
of functions fp^ s (x), which represent the positive and negative frequency 
solution at t — ► — oo. The field function can be decomposed as 

^) = E J ^^(^ + ^;k)a s (k)+^ + (x;k)6+(k)), (1) 

where we perform the sum over spin indices s. a s (bf) is the annihila- 
tion (creation) operator of the (anti)particle with the spin s and momen- 
tum k. If we demand that the creation and annihilation operators intro- 
duced in this way satisfy general anticommutational relations, then the 
Fock space can be constructed with the vacuum state defined at t — > — oo. 
However, following this procedure one can easily prove that Hamiltonian 
in the external time-dependent field will be non-diagonal (except the ini- 
tial state t — ► — oo). To diagonalize it one can use Bogoliubov time- 
dependent transformation to redefine the particle creation c s (t) and anni- 
hilation dj (t) operators with nontrivial equations of motion for them. The 
particle number in this formalism can be defined at every time moment as 
N ks (t) = fs(Kt) = (0|c+(k;i)c s (k;i)|0). 

1 The Bogoliubov transformations are named after N.N. Bogoliubov who introduced 
and applied these canonical transformations to determine integrals of motion for 
creation and annihilation operators [12] , 

2 As was noted in [5] the linear polarised field is not appropriate for a quantitative 
description of the laser pulse. However ,the linear polarised field was considered [2l[5] 
|6l [13] to obtain qualitative level results. 
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In this work, we will use another approach — the so-called "oscillator 
representation", which was earlier developed for bosonic systems \15\ 116]. 
The main idea is to construct the field function and to define the time- 
dependent creation and annihilation operators so that the Hamiltonian im- 
mediately becomes diagonal. Moreover it should take the form of the free 
field Hamiltonian at t — > — oo. The solution of equations of motion de- 
rived from Dirac equation for introduced operators will provide the wanted 
distribution function f s (k,t). 

This paper is organised as follows. In Section 2 the oscillator represen- 
tation formalism for fermions is given. In Section 3 the kinetic equation is 
derived. The applications of the model and particularly the fermion produc- 
tion in the early Universe are considered in Section 4. Section 5 summarises 
the article. Appendix contains the KE in the convenient form for numerical 
calculations. 



2. Oscillator (Holomorphic) Representation 

We start from the Lagrangian for electrodynamics 

L = 4>(iD^ - m(t))^, (2) 
where D a is the covariant derivative 

D a = d a + ieA a . (3) 

Let us assume that the function m(t) is an arbitrary smooth function 
of time. We also demand that the electric field is classical and space- 
homogeneous with a continuous A a (t). The classical approximation is sup- 
ported by the leading order of a large N approximation, see e.g. |10| . 

In order to derive the correct oscillator representation of the fermion 
field one can use the Foldy-Wouthuysen transformation |17j : 

HX) = I (^^ eiPXe ' 5(P,i) MP'*H(°) + bt(-P,t)v s (0)) , (4) 
^ {X) = /(2^ e " PX K(P>*) S '(°) +b.(-p,t)v a (0))e iS <P> t \ (5) 

where summation over spinor indices (s = 1,2) is implied, and spinors 
u(0),v(0) are defined as a solution of free Dirac equation in the particle rest 
frame: 



( 7o - l)u(0) = 0, 
(70 + lMO) = 0. 



(6) 
(7) 
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The unitary operator e lS ^ p '^ can be dependent both on time and mo- 
mentum. According to the scheme of oscillator representation e iS ( p ' t ' is 
defined so that the Hamiltonian density after the field decomposition 

H(p,t) = (a+(p,t)u s (0)+b s (-p,t)v s (0)) x 
e iS(p,t) (Pm + m)e iS(p,t) x 

(a s (p,t)u s (0) + b+(-p,t)v s (0)) (8) 

becomes diagonal, i.e. 

H(p,t) = (af(p,t)u s (0)+b s (-p,t)v s {Q))x 

u(P, t) x (a s (p, t)n s (0) + 6+(-p, t)v„(0)) , (9) 

with quasiparticle energy cv(P,t) = \/P 2 + m 2 defined by kinetic momen- 
tum P = p — eA. 

This condition can be rewritten as 

e «?(p,t)( P . 7 . + m ) e iS(p,t) = u (p,t), (10) 

with the following solution 

e iS(p,t) = ^(P,t)+m + P, 7 i (n) 
y/2u(P,t)(m + u;(P,t)) 

One can write down the explicit expression for S in a close form 



S(p,i) = i^, (12) 



where <ft is given by 



sin(2</>) = ^-J-, cos(2«/») = — — . (13) 
cj(P,tj cj(F,i) 

It is easy to convince that 

n(p^P) = e iS ^u(0), (14) 

u(-(p->P)) = e iS ( p '*)u(0), (15) 

where u(p),i;(p) are the solutions of free Dirac equation. Taking into ac- 
count these expressions (|14p . (|15|) we can rewrite the decomposition of Eqs. 
(JH), © in the form 

^) = / j^e^ {a s (p,t)u s (P) + bt(-p,t)v s (-P)} , (16) 
^) = /^ 7i e- i P x {a+(p,t)^(P)+6 s (-p,^(-P)}. (17) 
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Thus, the formal substitution p — * P in the free field decomposition results 
in ©, ©. 

The equations of motion for the creation and annihilation operators can 

be obtained from Dirac equations (dots denote the derivative with respect 
to time): 

<v(p,i) = -iooa s ,(p,t) - & u s , s a s (P,t) -E u s , s bf{-p,t), (18) 

6+(-p, t) = iub+(-p, t) - 9^6+(-p, t) - Z v s , s a s (p, t), (19) 
or in matrix form 

a(p,t) = -iua(p,t) - O u a(p,t) -E u b + (-p,t), (20) 

b+(-p, t) = iu;b + (-p, t) - @%+(-p, t) - E v a(p, t), (21) 
where Q u ' v , 3 U ' V are given by 

e^, s = u+(0)e- lS d e iS u s (0), (22) 

= vf,(0)e- iS d e iS v s (0), (23) 

S«, s = u+(0)e- iS d e iS v s (0), (24) 

^ = vU0)e- tS d e iS u s (0). (25) 
The straightforward calculation gives 

e u = e v = ie [PxE \ cr, (26) 
2u> (uj + m) 

S- = —E v = ^ + 7h Per - ^Ea, (27) 

2uj(uo + m) 2w ' 

where are Pauli matrices and E = — A is the field strength. 

In infinite past t — > — oo, the operators o s and b s should satisfy the 
classical anticommutation relations. However one can see that the second 
terms in the right part of the equations of motion (|20l [2T|) can modify the 
classical anticommutation relations in the presence of the time-dependent 
field. Indeed let us assume that 

K(k,t),a+(k',i)} =Il ss ,(k,t)5 3 (k-k'), (28) 

where the equation for unknown time-dependent matrix H ss ' is derived by 
using ([201 EED 

n(k,t) = [n(k,t),e u (k,t)]. (29) 

However, taking into account trivial symmetries of the matrix Tl ss i 

Il ss , = U s , s = IL* fs (30) 

and initial conditions in infinite past n ss /(k, t — > — oo) = 5 SS / one easily 
proves that n ss /(i) = const = 6 ss r, i.e. canonical commutation relations are 
valid even in time-dependent external field. 
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3. Distribution Function and Kinetic Equation 

In this section we will derive the equation of motion (or KE) for the 
one-particle correlator 

f s Ap,t) = (0|a+(p,t)a s (p,t)|0). (31) 

The distribution function is the diagonal components of f ss '(p, t) [9| [101, [18] . 
Further calculations will show that f ss i is coupled to one-particle correlator 
describing the antiparticle distribution g ss / and two anomalous correlators 

y ± r. 

g ss >(p,t) = (0|Zv(-p,i)6+(-p,i)|0), (32) 

V7s>M = (0|& s '(-P,*)a s (P,*)|0), (33) 

y+, (p,t) = (0\a+,(p,t)b+(-p,t)\0). (34) 

KE is obtained by differentiating f|31ti34|) with respect to time and subse- 
quent substituting the equations of motion (|2QII2ip and definitions (|31ti34j) . 
After some trivial transformation one gets: 

f = i[f,e]-y-E-Ey+, (35) 

y- = -2iuy~ +i[y-,e) + fZ-Eg, (36) 

y + = 2iu;y + + i[y + , 9] - gE + Ef, (37) 

g = i[g,@]+y-E + Ey + , (38) 

where = Q u /i = 0j<7j and E = E u = HjCXj are hermitian matrices. It is 
more convenient to introduce two hermitian matrices 

r + = \(y + + y-), (39) 
r" = l -(y + -y-) (40) 



and decompose \ = /> <?> r± hito U(2) basis (/ is a unitary matrix) 



X = XoI + XiVi, i = 1,2,3; (41) 
X0 = ^Tr X , Xi = ^Tr(x^). (42) 
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Finally in terms of (|41H42I) KE (|35H38|) takes the following form: 

f = -2Sr+, (43) 

f = -2[f x 0] + 2[r - x 3] - 2r+3, (44) 

90 = 2Hr+, (45) 

g = -2[g x 0] - 2[r" x S] + 2r+3, (46) 

r+ = 2ur Q + 3f - 3g, (47) 

r+ = 2o;r- - 2[r+ x 0] + 3(/ - 9o ), (48) 

r - = -2wr+ (49) 

r~ = -2ujt + - 2[r~ x 0] + [f x 3] + [g x 3]. (50) 

The KE (|43M50p allows some further simplifications if we assume zero 
initial conditions in infinite past t — > — oo for /q, f = (absence of particles), 

rj,r ± = (anomalous correlators turns to zero at in- vacuum), and go = 
l,g = (neutrality condition and absence of antiparticles). Summing (|4"3"]l 
and (|45p. (|44p and (|46p one finds the following integrals of motion: 

fo = 1-50, (51) 

f = -g. (52) 

Taking them into account KE equation (|43H50p is rewritten as 

f = -23r+, (53) 

f = -2[f x 0] + 2[r~ x 3] - 2r+3, (54) 

r+ = 2ujr +2ai, (55) 

r+ = 2u>r~ - 2[r+ x 0] + S(2/ - 1), (56) 

r = -2^r +, (57) 

r~ = -2wr + - 2[r - x 0]. (58) 

4. Applications 

4-1. Scalar-like interaction 
If the vector field vanishes Ai(t) = Q,Ei(t) = 0, we obtain 

= 0, (59) 



m 
2^ 



2 P (60) 



For this case the KE (I43M50P allows some simplification. First of all, the 
vector components r 1 * 1 , f, g will be collinear to p. Thus let us redefine 

r± = T^j^i (61) 
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and assuming zero initial conditions in the asymptotic past rewrite the KE 
in the form 

/o = 2Wv + , (62) 
t)+ = W(l-2/ )-2wt;-, (63) 
v~ = 2ujv + , (64) 

where W = ^V|p| and in the above equations the neutrality condition g = 
1 — /o was used. Applying zero initial conditions to /o and f we can conclude 
that the vector components of the distribution function for fermions and 
anti-fermions are identical zero and the evolution of /o component is given 
by §MM- 

So far we have not discussed any origin of the time dependence of the 
mass m(t). However any fermion interaction of the Yukawa type AL = 
g4>ipip in the mean field approximation for the scalar field (f) results in the 
time-dependent effective fermion mass m(t) = m + g{cp). Such a type of 
interactions is presented in the Higgs model, sigma model and can appear 
in some effective field theories (e.g. Walecka model |19j). 

Here we will concentrate on another important example where time- 
dependence of the effective fermion mass can be generated, namely, confor- 
mal field theory. 



4.1.1. Fermion production in the Early Universe 

For the sake of consistency let us recall some general and well-known 
facts of general relativity. 

The Lagrangian of the massive fermion field in general relativity is writ- 
ten as 

L = y/—g [iip(x)~f n (x)A n ifj(x) — m^j(x)ip(x)] , (65) 

where 

l n (*) = h n {a)1 a , (66) 
A„ = dn + hj^h^Yl^ (67) 

The rotation coefficients can be defined by means of tetrad as (for details, 
refer to 0[2O]) 

C^vp = {\nh^))h{y)rJtf^y (68) 

We will consider conformally static space-time. In this case the met- 
ric is conformally related to a static space-time. For a particular case of 
Robertson- Walker metric we have 

9»v = a 2 (r])hp U , (69) 
ds 2 = a 2 (r))(dri 2 - dx 2 ), (70) 
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Fig. 1. The time dependence of the number density for fcrmions with different 



masses 



2H, 



0.5, 1, 2 (from top to bottom). 



where a is the conformal factor, r] = J dt a(t) is the conformal time, 

is the Minkowski metric. Introducing new field variables tp c = a 1 / 2 (77)7/; we 

obtain 

L = a 2 (r]) [ii> c (x)^d^ c (x) - m(ri)ip c (x)'i() c (x)] , (71) 

where we define the effective fermion mass 771(77) = ma{rj). 

The cosmological creation of particles is considered in terms of the so- 
called conformal Universe [21] , where the volume of the Universe does not 
increase, while all masses, including the Planck mass, are scaled by the 
cosmic factor a(rj). We assume a stiff state scenario, where the mass 777,(77) 
is given by 

777(77) = "7/^/1 + 2^/77, (72) 

where Hi is initial data at the matter production instant. 

The kinetic equation transformed into convenient form for numerical 
calculations is given in Appendix. In Fig. 1 the evolution of the fermion 
number density (gj is the degeneracy factor) 

n = 2g f J /(P^) ( 73 ) 

for different ratios jjjj is shown. The final spectra of produced particles 
are far from exponential as seen from Fig. 2. The time evolution of specific 
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10" s| 1 

1 2 3 4 5 

q 

Fig. 2. The momentum dependence of the distribution function for fermions with 
different masses = 0.5, 1, 2 (from top to bottom) at the fixed moment of time 
r = 2Hit = 6; q = |q| is dimensionless momentum q = p/mj. 

entropy S/n is shown in Fig. 3, where entropy S is given by 

S = 2g f J ^3 [/(p,t)ln/(p,t) + (l-/(p,t))ln(l-/(p,t))]. (74) 

12. i 




1 2 3 4 5 6 

Fig. 3. The time dependence of the specific entropy S/n (|74p for fermions with 
different masses jjf- — 0.5, 1, 2 (from bottom to top). 
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4-2. Vector interaction 

In this section we consider some application of the model to the inter- 
action of the type AL = eA^tp^ij). 

4.2.1. Particle production in the external field with fixed direction 

Here we consider another simple example, particle production in the 
external time-dependent field with fixed direction. As in the previous sub- 
section, let us demand zero initial conditions for the particle fields in the 
asymptotic past. 

The vector potential is given by time-dependent amplitude and constant 
vector n: 

A = A(t)n. (75) 
The KE in this case also can be reduced^ to the system of three equations 

fa = 2Wv+, (76) 

v + = W{l-2f )-2uv-, (77) 

v~ = 2ujv + , (78) 

with W, defined by 



e|Ek/p1 + m 2 

w = |s 'r — ~w 2 — ' (79) 

v ± = =F(3r»f + ^-A ), (80) 
where the vector decomposition was used: 

P = An+i^pi, (81) 
I P-L I 

r ± = r ±n + -^-rf. (82) 
11 Ip±I 

The set of equations (|76M78p coincides with pioneering results obtained 
by Bialynicki-Birula, Gornicki, and Rafelski \22\ in the framework of the 
equal-time Wigner function. Also, for the one-component external field 
A = (0,0,^3(4)), the derivation of this KE in the framework of Bogoli- 
ubov quasiparticles can be found in [11]; some applications to particle 
production in strong laser fields \V'6\ [T4] as well as in the ultrarelativistic 
heavy- ion collision [H [23] in the Abelian dominance approximation [21] were 
also considered. 



3 We leave the proof to the reader. 
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5. Summary 

In the present paper, we use the oscillator representation [HJ [16] to 
derive the most general kinetic equation for fermions in time-dependent ex- 
ternal electro-magnetic field of arbitrary direction. Also particle production 
due to time dependence of fermion mass (as a result of Yukawa-like inter- 
actions in the meanfield approximation or conformally static space-time) is 
considered. 

The derived KE in the special cases (time-dependent mass, fixed direc- 
tion of the external field) coincides with KE obtained by time-dependent 
Bogoliubov transformation approach. However for the latter there is an 
unsolved problem of particle production in the vector fields of alterating 
direction. As shown in the present paper, the oscillator representation over- 
comes this difficulty and helps us to derive KE for fermions created in field 
of general polarization A = (Ai(t), ^(t), A^(t)). 

Some applications of derived KE are discussed. The detailed investiga- 
tion and numerical calculations are done for fermion creation in conformal 
cosmology. The solution of KE in this case has shown that the obtained 
momentum distribution function is far from equilibrium, therefore can be 
an argument in favour that the introduction of temperature is impossible 
without taking into account collision processes in created fermion gas. 

The consistent description of particle creation in the strong laser fields 
includes the effects of the field direction alteration. So far as the simplest 
model of the laser field only the linear polarised one is considered (e.g. 
[21 [5l O [13] ) . Our approach based on quantum field theory allows numerical 
calculation of particle production in the field of arbitrary polarisation. In 
this case the simplifications assumed in Section 4 are no longer applicable 
and one has to deal with the general system of equation (|53H58p . Preliminary 
numerical results on this subject are promising but out of the scope of the 
current paper and will be printed elsewhere. 
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Appendix 

We introduce dimensionless variables of the time r, momentum q and 
the parameter 7/ according to the following formulas: 

r = 2H IV , (83) 
q = p/mj, (84) 

•» ^ ^ < 85 > 

The kinetic equation (I62H64P in this variables takes the form 

^/o = 2Wv+, (86) 

^- v + = W(l - 2/o) - 2Cj 1iV ~, (87) 
ar 

d , 



where 



* = i (q2 + 1 | q |)VTTT - (89) 

w = Vq 2 + 1 + T. (90) 
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